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Abstract 

We discuss the extension of the oscillator-basis J-matrix formalism on the case of 
true A-body scattering. The formalism is applied to loosely-bound n Li and 6 He 
nuclei within three-body cluster models 9 Li + n + n and a + n + n. The J-matrix 
formalism is used not only for the calculation of the three-body continuum spectrum 
wave functions but also for the calculation of the S"-matrix poles associated with 
the n Li and 6 He ground states to improve the description of the binding energies 
and ground state properties. The effect of the phase equivalent transformation of 
the n—a interaction on the properties of the 6 He nucleus is examined. 



1 Introduction 



In this paper we discuss a unified J-matrix approach to many-body systems. 
The approach links the quantum scattering theory formalism with traditional 
variational methods of nuclear theory based on the wave function expansion in 
the harmonic oscillator function series. The formalism extends the variational 
methods on the case of the continuum spectrum states. On the other hand, the 
variational description of the discrete spectrum states is improved by the use 
of the methods of the scattering theory. The J-matrix approach is supposed to 
be a valuable extension of the shell model and other nuclear structure models. 
Below we examine in detail the application of the J-matrix approach to the 
light weakly bound nuclei n Li and 6 He within the three-body cluster model. 
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Weakly bound nuclei are extensively studied now both theoretically and exper- 
imentally. New methods are developed in the nuclear structure theory provid- 
ing accurate description of the charge and matter distribution in the nuclear 
halo, allowing for the coupling to the continuum spectrum states, etc. At the 
same time, the traditional language of nuclear theory, i. e. the nuclear shell 
model, is based on the oscillator basis expansion of the wave functions and 
seems to be inadequate for the the continuum or weakly bound states with 
the wave functions slowly decreasing at large distances since the oscillator 
functions are rapidly decreasing asymptotically. We however believe that the 
traditional oscillator-basis nuclear theory can be successfully developed to face 
the above challenges by implementing the J-matrix formalism. 

Within the J-matrix formalism, the continuum spectrum wave function is ex- 
panded in infinite series of L 2 functions. The J-matrix method was initially 
proposed in atomic physics [1,2] and shown to be one of the most efficient 
and precise methods in calculations of photoionization [3,4,5] and electron 
scattering by atoms [6] . In nuclear physics the same approach has been devel- 
oped independently [7,8] as the method of harmonic oscillator representation 
of scattering theory. This method has been successfully used in various nuclear 
applications allowing for two-body continuum channels, e. g. nucleus-nucleus 
scattering has been studied in the algebraic version of RGM based on the J- 
matrix formalism (see the review papers [9,10]); the effect of A and neutron 
decay channels in hypernuclei production reactions has been investigated in 
Refs. [11,12], etc. 

In this paper we study exotic neutron-excess nuclei 11 Li and 6 He in the three- 
body cluster models n Li = 9 Li + n + n and 6 He = a + n + n. The two-neutron 
separation energy in both of these nuclei is small compared with the neutron 
separation energy in the 9 Li and 4 He clusters; as a result the wave function 
of the pair of neutrons decreases slowly with distance and the rms radius of 
the two-neutron distribution is large compared with the rms radius of the 9 Li 
or 4 He core (the so-called two-neutron halo). We note that n Li and 6 He are 
the so-called Borromean nuclei, i. e. none of the two-body subsystems in the 
cluster decompositions 11 Li = 9 Li + n + n and 6 He = a + n + n has a bound 
state. We employ the three-body extension of the J-matrix formalism for the 
description of the three-body decay channels and develop this approach to 
improve the description of the ground state properties. 

Three-body decay channels are much more informative than the two-body 
ones, however the complete description of a three-body decay is a very com- 
plicated problem. Nevertheless in some cases, one can use the so-called demo- 
cratic decay approximation for the successful analysis of the experimental data 
or for obtaining reliable theoretical predictions. 

Generally, a three-body continuum spectrum wave function in the asymptotic 
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region is a superposition of the components describing two- and three-body 
decay channels [13]. The democratic decay approximation accounts for the 
three-body channel only; two-body channels associated with the appearance 
of bound two-body subsystems and 1 non- democratic 1 subdivision of the system 
of the type (2+1), are not allowed for in the approximation. Therefore, this 
approximation is valid for the study of a three-body system only in the case 
when all two-body decay channels are closed and the only open channel is a 
three-body one; in other words, this is the case when none of the two-body 
subsystems has a bound state. Hence the democratic decay approximation is 
adequate for the study of Borromean nuclei 6 He and n Li at least at small 
enough excitation energies (less than the single nucleon binding energy in the 
9 Li and 4 He clusters). 

Generally speaking, if the three-body channel is the only open one, than the 
system can have two types of asymptotics [13]. One of the asymptotics cor- 
responds to the situation when one of the particles is scattered by another 
and the third particle is a spectator. From the general physical point of view, 
this asymptotics is supposed to be of little importance for a nuclear system 
excited in some nuclear reaction and decaying via a three-body channel. The 
alternative type of the asymptotics is a superposition of ingoing and outgoing 
six-dimensional spherical waves, it corresponds to the situation when the de- 
caying system emits (or/and absorbs) three particles from some point in space. 
Allowing for the three-body asymptotics of this type only is a quintessence 
of the democratic decay approximation; sometimes this approximation is also 
referred to as true three-body scattering or 3 ^ 3 scattering. 

Nuclear structure studies in the framework of the democratic decay approx- 
imation have been started in 1970th by R. I. Jibuti and collaborators. The 
review of the results obtained by the Tbilisi group can be found in Refs. [14,15]. 
Recently a considerable progress has been made in both theoretical and ex- 
perimental studies of democratic decays (see reviews in Refs. [16,17,18,19]). 

The generalization of the J-matrix formalism on the case of 3 — > 3 scattering 
(and on a more general case of N — > iV scattering) was suggested in Ref. [20]. 
First it was successfully applied in the study of monopole excitations in 12 C 
nucleus in the cluster model 12 C = a + a + a in Ref. [21]. Later by means of 
this approach we studied AA hypernuclei [22]. Our first attempts to employ 
this approach in the studies of n Li and 6 He nuclei in the three-body cluster 
models can be found in Refs. [23,24,25,26]. 

We use the J-matrix 3^3 scattering formalism for the construction of the 
three-body continuum spectrum wave functions of n Li and 6 He nuclei excited 
in nuclear reactions. However, probably more interesting is the application of 
this formalism to the study of the ground states of these nuclei treated as 
three-body systems 9 Li + n + n and a + n + n. The idea is the following. 
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By the J- matrix formalism we calculate the 3 — > 3 scattering S- matrix. This 
calculation may be extended on the complex momentum plane. Hence we can 
locate numerically the S-matrix poles associated with the bound states. The 
S-matrix pole calculations improve the variational results for binding energies 
obtained by the pure diagonalization of the truncated Hamiltonian matrix. 
Knowing the binding energies, we calculate the bound state wave functions by 
means of the J-matrix formalism. We suppose that the approach based on the 
S-matrix pole calculation will be useful not only within the cluster models; 
the implementation of the A — > A scattering S-matrix calculation within the 
J-matrix formalism is believed to improve essentially the traditional nuclear 
shell model. 



The first calculation of the S-matrix poles within the J-matrix formalism was 
performed in Ref. [27] where resonances in a— a scattering were investigated. 
The S-matrix poles in He atom and H~ ion were calculated in Ref. [4,5] in the 
J-matrix model suggested by Broad and Reinhardt [3]. The resonance energies 
and widths were reproduced with a high accuracy. The S'-matrix poles asso- 
ciated with the bound states were also calculated in Ref. [4,5]. The binding 
energies were reproduced with a very high accuracy; we note also that a very 
large number of bound states was obtained by means of the S'-matrix pole 
calculations: it was not only much larger than the number of bound states 
obtained variationally (by the diagonalization of the truncated Hamiltonian 
matrix), but even much larger than the rank of the truncated Hamiltonian 
matrix. This approach not only provides an essential improvement in calcu- 
lations of binding energies, but also makes it possible to calculate the ground 
state wave function with the correct asymptotics. It is very important for halo 
nuclei like n Li and 6 He due to the slow decrease of the wave function in the 
asymptotic region that cannot be reproduced by a superposition of a finite 
number of oscillator functions obtained in variational calculation. As a result, 
the rms radius, electromagnetic transition probabilities and other observables 
are improved essentially. 



We also discuss a phase-equivalent transformation suggested in Ref. [26] that 
is used to obtain families of phase-equivalent n— a potentials. The effect of the 
phase-equivalent transformation of the n— a interaction on the 6 He properties 
is examined. 



We start the discussion with a brief sketch of the N — > N scattering J-matrix 
formalism suggested in Ref. [20]. 
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2 J-matrix formalism with hyperspherical oscillator basis 

The wave function \P of a system of A particles is generally a function of 
A coordinates of individual particles ty The center of mass motion can be 
separated and eliminated; as a result we can treat the wave function \l/ as 
a function of A — 1 Jacobi coordinates Within the democratic decay ap- 
proximation, it is natural to employ the hyperspherical harmonics formalism 
(see, e. g., [15,28]). The 3A — 3 independent variables £j are rearranged in this 

M-i 2 

formalism: the so-called hyperradius p = W ^ is introduced and all the 

rest — 4 variables are the angles on the (3 A — 3)-dimensional sphere. The 
wave function is searched for in the form 

^ = ]Tdr^r, (1) 
r 

where 

= V Kl = $^ 7 (p) 0fc 7 (ft), (2) 

f2 is the set of angles on the (3 A — 3)-dimensional sphere, K is the so-called 
hypermomentum [physically K is the angular momentum in the (3 A — 3)- 
dimensional space; K > L where L is the orbital angular momentum] and 7 
stands for all the rest quantum numbers labelling the hyperspherical function 
£2^ 7 (f2), the multi-index T = {K,^}. Various analytic expressions may be 
found for the hyperspherical functions £%- 7 (f2) in textbooks (see, e. g., [15]). 

In the two-body case, p = r = |t*i — T2\ is the distance between the particles 
and $k 7 (f2) becomes the usual spherical function Ylm(&) with K and 7 play- 
ing the role of the angular momentum L and its projection M, respectively. 
Therefore equation (2) is a generalization on the A-body (A > 3) case of the 
conventional wave function used in the study of spherically symmetric two- 
body systems. However, contrary to the case of two-body systems with central 
interaction which impose the spherical symmetry, the hypermomentum K is 
not an integral of motion: the two-body interactions in the A-body system 
couples the states with different values of the hypermomentum K and we 
have a set of coupled equations for the hyperradial functions $x 7 (p). This set 
of equations has the same structure as the one describing two-body systems 
with non-central forces; from the point of view of scattering theory, this set 
of coupled equations is formally equivalent to the one describing multichannel 
scattering in the system with the hyperspherical channels Y. 

We introduce the hyperspherical oscillator basis 

|«r) = \ K Ki) = m kK { p ) 0fc 7 (n), (3) 
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where the hyperradial oscillator function 

«w( P ) = <s% u>) = P -« M - 4 >/ 2 t KK ( P ), (4) 

is the associated Laguerre polynomial and 

In the two-body case (A = 2), «5f is equivalent to the orbital angular momen- 
tum L, J2? = L, and Eqs. (3)-(6) define the conventional oscillator basis. In the 
many-body case (A > 3), Eqs. (3)-(6) define the A-body oscillator functions, 
i. e. the eigenf unctions of the A-body Schrodinger equation with the potential 
energy 

U=\ E H) 2 = f P\ (7) 

where m« is the mass of the ith particle, R is the center-of-mass coordinate 
and ui is the parameter of the oscillator basis. The corresponding eigenenergy 
is 

E KK =(n+ tkj, (8) 



2 

where the number of oscillator quanta 

3A — 6 

N = 2k + K = 2k + ^ — . (9) 

The hyperspherical oscillator basis (3) is orthonormalized: 

(nK^n'K'i) = <W<5rr', (10) 



or 



/ <kWk'k(p) dp = <W- (11) 



The wave function \1> is expanded in the hyperspherical oscillator function 
series, 

m = ]T (kK^) IkK'j), (12) 

K,K,J 

and the Schrodinger equation takes the form 

E (kK^H - E\K'K'i) {K'K'i\^) = 0, (13) 

k',K',Y 
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where E is the energy, the Hamiltonian H = T + V, the potential energy 
V is usually a superposition of two-body interactions V^-, V = Vij, an d 

i<j 

T is the A-bodj kinetic energy operator. Generally, the Hamiltonian matrix 
{kK^\H\k! 'K' '7') is infinite. However, within the J-matrix formalism we trun- 
cate the potential energy matrix; it is most natural to define the truncation 
boundary through the number of oscillator quanta, i. e. we suppose that 

{kK^\V\k' K'^') — if 2k + K>N or 2k + K' > N . (14) 

The kinetic energy matrix (K,K / y\T\K > K' / -f') is tridiagonal, 



{KK 1 \T\K>K'i) = f ^-5 KK ,5 11 , 



\2n + J?+-)5 KK , 



- y (n + 1) (n + J&? + £) <W,*' 



• (15) 



The A-body hyperspherical J-matrix formalism is very close to the conven- 
tional J-matrix formalism with oscillator basis [2,29] used in multichannel 
two-body problems. The main difference is that Jzf entering Eq. (15) is half- 
integer if the number of particles A is odd [see (6) and note that K is always 
integer] , in particular, Jzf is half-integer in the three-body case discussed here. 
The general J-matrix oscillator-basis solutions of the free Schrodinger equa- 
tion 

{kK-i\T - E\K'K'i) {n'K'i\m) = (16) 

k',K',Y 

that may be used in the case of arbitrary Jzf, were suggested in Ref. [20]. The 
regular solution is 

^W^ r( , i + ^ + 3/ 2 ) ^-'' /2 ^V). (17) 
the irregular solutions are 

C KK (q) 



■lq y p J S 0K (q') S K K(q') ^ 



nS 0K (q) 

2q 7S 0K (q')S KK (q' 



q 2 — q' 2 



I 



KS OK (q)J o q 2 -q' 2 + t0 



dq', 



(18) 



(19) 



2g f S 0K (q') S KK (q') d ^ 



*S OK (q)J o q 2 -q' 2 -t0 



(20) 
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2E 

where q = \ — , V. P. in Eq. (18) indicates the principal value integral, 
V nw 

+iO and — iO in Eqs. (19)-(20) show how the poles of the integrand should be 
treated. The solutions (17)-(20) are not independent, the following relation 
between them is valid: 

C { ^(q) = C KK {q) ± iS KK {q). (21) 

However, any pair of the solutions (17)-(20) can be used to construct an 
arbitrary solution. 

The regular solution (17) is just the hyperradial momentum-space oscillator 
function. The irregular solutions (18)-(20) are more complicated. They were 
analyzed in detail in Ref. [20]. In the general case, they can be expressed 
through Tricomi function. However in the case of even A, the irregular solu- 
tions can be simplified and expressed through the confluent hypergeometric 
function. The physical meaning of the solutions (17)-(20) is clear from the 
following expressions 



00 2 ( nJ£\ 

Y, S KK(q)*K K (p) = v^^+iM^V- sin Up - — I , (22) 



K=0 



E ^(?) ^ - # (gp) ^ l| - coslgp-^-J, (23) 



K=0 

OO 



E^WWrf-^^M-^ e<--^), (24) 



K=0 

t^^MA^-iVTpH^p)^ 2 - e-<--*). (25) 

Here J a (x), N a (x) and H^' 2 \x) are Bessel, Neumann and Hankel functions, 
respectively. 

In the case of continuum spectrum (E > 0), the oscillator representation wave 
function (kK^\^^ t )) = (fip|\l/( r )) in the channel T is of the form 

W^) = \ (5 Tr C$(q) - c£>(q) [S]nv) (26) 

in the 'external region' k > k t , where n r = (N — K)/2 is the potential energy 
truncation boundary in the channel P It is supposed that the ingoing spherical 
wave is present in the channel P only while the outgoing spherical waves are 
present in all channels; [S]rr ; entering Eq. (26) is the matrix element of the 
hyperspherical S'-matrix. The S'-matrix [S] can be calculated by the following 
formula [20]: 

[S] = [A s ] _1 [B s ] , (27) 
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where the matrix elements of the matrices [A s ] and [B s ] are 

[A% r = (« r ,r'|^ r + 1, r> C^ K (q) - 5 TVI C^ K (q), (28) 

[B*] rT = (« r ,r'|^ r + 1, r> c^ K ( q ) - 5 rr c^M; (29) 

the matrix elements of the matrix \&\ proportional to the discrete analog of 
the P-matrix (see Ref. [29] for details) are 

{kT\@>\k v , + = ^\%k t ,V) {k v ,V\T\k t , + 1,V), (30) 

the kinetic energy matrix elements (/t r T'|T|/t r / + 1, r") are given by Eq. (15), 
and the matrix [^3] = [H — E] 1 defined in the truncated model space spanned 
by oscillator functions (3) with k < Kr — (N — K)/2 in each channel T, has 
the matrix elements 

(«rW*) = E™ (3D 

A ^ 

In Eq. (31), E\ are eigenenergies and (fiT|A) are the respective eigenvectors 
of the truncated Hamiltonian matrix, i. e. E\ and (kT|A) can be found by 
solving Eq. (13) supposing that k < kt and k' < kt 1 - The results obtained 
by diagonalization of the truncated Hamiltonian matrix we shall refer to as 
variational results. 



After performing variational calculation, we obtain the S-matrix with the help 
of Eqs. (27)-(31) and the oscillator representation wave function (fiT|\I/ <T )) in 
the external region k > k t by the formula (26). The oscillator representation 
wave function (kT\^^) in the 'internal region 1 k < k t can be now calculated 

as 

(«r|*( r ')) = J2 (Kr|^|Krv, + i,r") (« r » + i,r"|* (r ' ) ), k<k t . (32) 

r" 



Equation (27) can be used in the complex momentum plane; the bound state 
energies are associated with the S'-matrix poles that can be found by the 
numerical solution of the obvious equation 

det [A s ] = 0. (33) 

The matrix [A s ] in Eq. (33) is the extension on the complex momentum plane 
of (28); for the the bound states with E = \q 2 hw < its matrix elements are 

[A*] rT = (k ti T'\0>\k t + 1, r> C^ K {q) - 5 TT , C^ K (q), (34) 

where 

C&(q) = C^(i\q\)- (35) 
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The relation (26) for the oscillator representation wave function in the external 
region k > K r should be replaced for bound states by 

(«r|*> = -[S] rr cS(2); (36) 

the multipliers [S]rr are obtained by the numerical solution of the equation 

[S]rr<k?) = E («rr|^|«iv + l,r') [SfrrCQ^q) (37) 
r' 

[we note that Eq. (33) is the condition of solvability of Eq. (37) for the bound 
states]. Equation (32) can be used for the calculation of the bound state os- 
cillator representation wave function (fiT|\l/ < ' r )) in the internal region k < n r . 
Since the set of [S]rr can be obtained from Eq. (37) up to a common multiplier 
only, the bound state wave function should be normalized numerically. 

The bound state energies and wave functions obtained by numerical calcula- 
tion of the S'-matrix poles, will be refered to as J -matrix results. 

In practical applications, it is useless to allow for all possible channels T = 
{K, 7} in the external region N = 2k + K > N. We start the calculations 
allowing for the channels T with few minimal possible values of K in the 
external region only, i. e. we allow for the channels with all possible values 
of K < K tr and all possible 7 values for a given K in the external region. 
Therefore the summation over V" in (32) is restricted to these allowed channels. 
However all possible channels with K < N are allowed for in the internal 
region N = 2k + K < N . As a next step, we increase the value of K tr and 
allow for more channels T in the external region. The convergence of all results 
(binding energies and other bound state observables, transition probabilities, 
etc.) with K tT is carefully examined. The convergence is usually achieved at 
small enough values of K tT (much less than N). Such converged with respect 
to K tI results are discussed below. 



3 Application to xl Li nucleus 

The n Li nucleus is studied in the three-body cluster model 9 Li + n + n. The 
three-body cluster model for n Li was first suggested in Ref. [30]. The n Li 
three-body binding energy or equivalently the n Li two-neutron separation 
energy is 0.295 ± 0.035 MeV [31] that is much less than the excitation energy 
of the lowest excited state in the cluster 9 Li. As it was already noted, n Li is 
a Borromean nuclei, i. e. none of the two-body subsystems 9 Li + n and n + n 
has a bound state. 

The two-body potentials are needed for the investigation of the three-body 
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system n Li = 9 Li + n + n. Unfortunately the information about the n— 9 Li 
interaction is scarce. As a result, in Ref. [30] the n— 9 Li potential was phe- 
nomeno logically parametrized and fitted to the 11 Li binding energy. Therefore 
it is reasonable to use a simplified three-body model with simplified interac- 
tions to avoid computational difficulties; a microscopic extension of the model 
seems to be useless due to uncertainties of the n— 9 Li interaction. 

We employ the n—n and n— 9 Li potentials of Ref. [30]. The Gaussian n— n 
potential 

U(r) = -V e- (r/i?)2 (38) 

was fitted to the low-energy singlet (S — 0) s wave phase shifts, its depth 
V = 31 MeV and its width R — 1.8 fm. The n— 9 Li potential is of a two- 
Gaussian form, 

U(r) = -V x e- {r ' R ^ - V 2 , (39) 

with V 1 = 7 MeV, R ± = 2.4 fm, V 2 = 1 MeV, and R 2 = 3.0 fm. It is sup- 
posed that two valent neutrons in 11 Li are in a singlet spin state. The total 
angular momentum J in 11 Li results from the coupling of the spin | of the 9 Li 
cluster with the relative motion orbital angular momentum L; however we do 
not make use of the spin-flip operators in our investigation, therefore we can 
disregard the spin variables and exclude them from the wave function. The 
11 Li ground state orbital angular momentum L is supposed to be zero. We can 
define the channels T as T = {K, l x , l y , L} where l x and l y are orbital angular 
momenta corresponding to the Jacobi coordinates 



- = ^-^_(r.-r 1 ), (40) 



_ juj m 3 (m 1 + m 2 ) ( myTi + m 2 r 2 \ 
V h mi + m 2 + m 3 \ mi + m 2 J ' 

the angular momenta l x and l y are coupled to the orbital angular momentum 
L. The wave function is given by the general formula (12). The oscillator basis 
(3) in our case takes the form 

|«r) = \nKlJ y : LM) ee (p\ K Kl x l v : LM) = M K , K (p) & K j(ty, (42) 

where the hyperspherical functions 

& K ^n) ee & KlxlyLM (n) 

= N l £ ly cos'"a sin'»a P^ +1/2 '^ +1/2 (cos 2a) 

Tn x lyTfly \LM)Y lxtmx (x)Y lvtmy (y), (43) 

m x +m y =M 

a = tan _1 (y/a;), Yi m is the spherical function, {l x m x l y m y \LM) is the Clebsch- 
Gordan coefficient, P"'P(x) is Jacobi polynomial [32], and the normalization 
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factor 

2n\ (K + 2) (n + l x + l y )\ 
T(n + l x + 3/2)T(n + l y + 3/2) ' 

The basis (42) is useful for the calculation of three-body decays. However it is 
not convenient in the calculation of the matrix elements of two-body potentials 
Vy. These matrix elements can be easily calculated in the basis 

\n x l x n y l y : LM) = (xy\n x l x n y l y : LM) 

= {l x m x l y m v \LM) (x\n x l x m x ) (y\n y l y m y ), (44) 

m x +m y =M 

where (x\n x l x m x ) and (y\n y lym y ) are the conventional three-dimensional os- 
cillator functions depending on the Jacobi coordinates x and y, respectively. 
The unitary transformations relating the basises (42) and (44), 

\kKI x I v : LM) = ^ {n x l x n y l y \hiK) \n x l x n y l y : LM), (45) 

^x i^y 

\n x l x n y l y : LM) = ^ (n x l x n y l y \KK) \nKl x l y : LM), (46) 

K,K 

are discussed in detail in Refs. [15,28]. The transformation coefficients 
(n x l x n y l y \K,K) are non-zero if only the the basis functions (42) and (44) are 
characterized by the same number of oscillator quanta N, i. e. when 
2n x + l x + 2n y + l y = 2k + K. 

We suppose that the 1st and the 2nd particles are neutrons and the 3rd particle 
is the 9 Li cluster (or the 4 He cluster in the next section). Therefore the Jacobi 
coordinate x is proportional to the distance between the valent neutrons and 
the Jacobi coordinate y is proportional to the distance between the cluster 
and the center of mass of two valent neutrons. In this case the basis (44) can 
be used directly for the calculation of the n— n interaction matrix elements. 
For the calculation of the matrix elements of the n— cluster interaction, we 
renumerate the particles assigning number 1 to the cluster and numbers 2 and 
3 to the valent neutrons, and use Eqs. (40)-(41) to define another set of Jacobi 
coordinates x' and y'. The n— cluster interaction matrix elements can be di- 
rectly calculated now in the basis (44). The unitary transformation of the wave 
functions associated with the switching from one set of Jacobi coordinates (x 
and y) to another (x' and y'), is discussed in detail in Refs. [15,28]. 

The n Li ground state energy dependence on the oscillator basis parameter 
hoo is shown on Fig. 1 for different values for the truncation boundary N. 
The dependence is typical for variational calculations and has a minimum 
at hu PS 6.5 MeV (the best convergence hu value). The hio ground state 
energy dependence obtained in variational and J-matrix calculations are of 
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Fig. 1. The 11 Li ground state energy vs the oscillator basis parameter huj in the 
variational (left panel) and J-matrix (right panel) calculations for different values 
of the truncation boundary N. The horizontal dashed line depicts the convergence 
limit for the ground state energy. 

the same type. However this dependence is much less pronounced in the J- 
matrix calculations (note a different scale on the left and right panels of Fig. 1). 
Therefore the J-matrix calculation is much less sensitive to the choice of the hu> 
value. The J-matrix ground state energy results are better than the variational 
ones for any value of hcu and any value of N. 

This is clearly seen on Fig. 2 where we present the convergence with N of the 
11 Li ground state energy and rms matter radius 



(r 2 ) l/2 = ^j^(r 2 c ) + (p 2 ), (47) 

where the n Li mass number A = 11 and (r^) 1//2 is the rms radius of the 9 Li 
cluster. The variational ground state energy decreases monotonically due to 
the variational principal; the variational rms radius consequently monotoni- 
cally increases. We note here that the variational principle is not applicable to 
the J-matrix calculations where the infinite oscillator basis is allowed for and 
the results obtained with different N values differ by the rank of the potential 
energy matrix. The J-matrix results demonstrate a much more interesting 
dependence on N than the variational ones. For the hcu values close to the 
best convergence value of fko — 6.5 MeV, there is a staggering of the J-matrix 
binding energy and rms radius with N (Fig. 2, left panel). The amplitude 
of the staggering decreases with N and the results converge rapidly. It is in- 
teresting that the binding energy and the rms radius converge to the values 
lying between the results obtained with two subsequent values of N . Hence we 
obtain not only the lower bound binding energy and rms radius estimates as 
in variational calculations but both lower bound and upper bound estimates 
for these observables. These estimates are presented in Table 1 together with 
experimental data and the results of calculations of other authors. 
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Fig. 2. Convergence with TV of the n Li three-body ground state energy (a) and 
rms matter radius (b) for hu>=6.5 MeV (left panel) and Huj=20 MeV (right panel). 
Variational and J-matrix results are shown by triangles up and triangles down, 
respectively. The convergence limits are shown by the dashed line. 



Table 1 

The n Li two-neutron separation energy E(2n) and rms matter radius (r 2 ) 1 / 2 ob- 
tained in the calculations with tvjj = 6.5 MeV together with the results of theoretical 
studies of Ref. [33] and experimental data. 



Approximation 


E{2n), MeV 
N = 38 N = 40 


(r 2 ) 1 / 2 , fm 
N = 38 N = 40 


variational 
J-matrix 
Ref. [33] 


0.326 0.327 
0.335 0.336 
0.3 


3.176 3.189 
3.336 3.348 
3.32 


experimental 
data 


0.247±0.080, Ref. [34] 
0.295±0.035, Ref. [31] 


3.10±0.17, Ref. [35] 
3.53±0.10, Ref. [36] 
3.55±0.10, Ref. [37] 
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Fig. 3. The dominant K = and X = 2, Z x = = compo- 

nents of the 11 Li ground state wave function in the oscillator representation, 
(k, K = 0, l x = 0, l y = : 00 | *) (circles) and (k, K = 2, l x = 0, l y = : 00 | *) 
(squares), obtained in the calculations with N = 40 and Hco = 6.5 MeV. The bold 
and empty symbols are the variational and the J-matrix results, respectively. 



If t\oj differs much from the best convergence value, we have the conventional 
variational-type monotonic convergence of the J-matrix ground state energy 
and rms radius (Fig. 2, right panel). The results converge slower with N in 
this case, however much faster than in the variational calculation. 

The structure of the ground state wave function is the following. The total 
weight of the l x — l y — components is 99.15%, in particular the K = 
component contribution is 94.4% and the K — 2, l x — l y — component 
contribution is 3.3%. This result agrees well with the results of other authors 
who used the three-body cluster model of n Li and the two-body interactions 
of the same type. In particular, the contribution of the l x — l y — components 
was estimated in Ref. [30] as 99%. According to Ref. [33], the weight of the 
dominant K — 0, l x — l y — wave function component is 95.3%, and the 
contribution of this component together with the next K — 2, l x — l y — 
component is 98.4%. The dominant components of the n Li ground state wave 
function in the oscillator representation, (k, K — 0, l x — 0, l y — : 00 | ^) and 
(«, K — 2, l x — 0, l y — : 00 | \&), are shown in Fig. 3. It is seen from the figure 
that the \ow-k dominating oscillator components are well reproduced in the 
variational calculation. However as k approaches the truncation boundary Kr, 
the variational calculation underestimates (k, K = 0, l x = 0, l y = : 00 | ^F) 
and (k, K = 2, l x = 0, l y = : 00 | essentially. 

A very interesting problem is the problem of low-energy electromagnetic tran- 
sitions in n Li and other neutron-rich loosely-bound nuclei. It is supposed that 
the El transition strength is strongly enhanced in such nuclei at small enough 
energies (the so-called soft dipole mode). We calculate the reduced El transi- 
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tion probability 



dB(El) 
dE 



2J. + 1 



1 



Y,\(Jf\\^(El)\\Ji) 



(48) 



where 




Zen / h 2m + M c 



A V u 2mM c 



(49) 



y is the Jacobi coordinate proportional to the distance between the 9 Li cluster 
and the center of mass of two valent neutrons, m is the nucleon mass and M c 
is the 9 Li cluster mass, e is the electron charge, the number of protons in n Li 
Z = 3, the number of valent neutrons n = 2, Jj = 3/2 and J/ are the total 
angular momenta of the ground (initial) and final states, respectively. We note 
here that within our three-body cluster model of n Li we calculate only the 
so-called cluster El transition strength associated with the relative motion 
of two neutrons and the 9 Li cluster; the total El transition strength includes 
additionally excitations of nucleons forming the 9 Li cluster that manifest itself 
at higher energies. Since at low energies the only open decay channel is the 
three-body one, n Li -> 9 Li + n + n, the democratic decay approximation is 
well justified for the calculation of the final state wave function. 



to compare the results obtained in the following approximations: 

VV — the variational calculation of both the ground and final states; 
VJ — the variational calculation of the ground state and the J-matrix cal- 
culation of the final three-body continuum state; 
and 

JJ — the J-matrix calculation of both the ground and final three-body 
continuum states. 

In the VV approximation, we obtain a discrete spectrum of excited states and 
cannot calculate the final state wave function at an arbitrary given positive 
energy. As a result, the El transition strength differs from zero at the energies 

belonging to the discrete spectrum of final states and — — — - has a form of a 

a E 

number of (^-function peaks shown in the figure by vertical lines with a cross 
at the end. The first peak at the excitation energy of approximately 2 MeV 
exhausts nearly 80% of the cluster energy-weighted sum rule [38] 



where E is the ground state energy. This is clearly seen from Fig. 5 where we 



The results of the 



dB(El) 
dE 



calculations are presented in Fig. 4. It is interesting 




(50) 
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Fig. 4. Reduced El transition probability 
dt3 dE^ m Vertical lines with cross at 
the end, dashed and solid lines were ob- 
tained in the VV, VJ and JJ approxima- 
tions, respectively, with hu> = 6.5 MeV 
and N = 20 for the ground state and 
N = 21 for the final state. Dash-dot and 
dash-dot-dot lines are the calculations of 
Ref. [39] and the parameterizations of ex- 
perimental data of Ref. [40], respectively. 

present the plot of the function 



Fig. 5. Cluster El energy-weighted sum 
rule in n Li. Dash-dot, dashed and solid 
lines were obtained in the VV, VJ and 
JJ approximations, respectively, with 
Huj = 6.5 MeV and N = 20 for the 
ground state and N = 21 for the final 
state. 



Cclust 
°E1 



(51) 



The 2 MeV peak should be related with the soft dipole mode in n Li in the 
VV approach. The energy of the soft dipole mode obtained in the VV ap- 
proximation differs essentially from the one of the experimental El strength 
low-energy maximum. 

In the VJ approximation, the final energy wave function can be calculated at 

, r , c , dB(El) 

any given positive energy. Instead of sharp o-peaks, we have a smooth — 

d E 

curve with a maximum shifted to a lower energy that is closer to the energy of 
the experimental El strength maximum. This is clearly the result of allowing 
for the continuum spectrum effects in the final state wave function. However 
the external 'asymptotic' part of the model space with 2k + K > N is not 
allowed for in the ground state wave function and hence the final state wave 

dB(El) 



function components with 2k + K > N + 2 do not contribute to 



dE 



The contributions to 



dB(El) 
dE 



of the external asymptotic part of the model 



space of both the ground state and the final state wave functions, are com- 
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Fig. 6. The dB }^ convergence with N at the energy E = 0.5 MeV where all 
components with the number of oscillator quanta 2k + K < N are allowed for in the 
JJ calculations with fojj = 6.5 MeV and N = 20 for the ground state and TV = 21 
for the final state. 

pletely accounted for in the JJ approximation. These contributions are seen 
from Figs. 4 and 5 to be of great importance. They shift the 



dE 



maximum 



to lower energy and change the shape of the El strength function. We note 
here that due to the slow decrease in the asymptotic region of the wave func- 
tion of the loosely-bound state, it is needed to allow for a very large number 
of components in the asymptotic part of the model space in the El strength 

dB(El) 



calculations. To calculate 



dE 



in the low-energy region with high accu- 



racy, we should allow for all components with the number of oscillator quanta 

2k + K < N where iV is of the order of 1000. The convergence of — j - 

d E 

with iV at the energy E = 0.5 MeV, is shown in Fig. 6. We note here that the 
Nhu oscillator function has the classical turning point at 



h 



l2Nh 



f52) 



Therefore supposing, say, iV = 500, we allow for the distances up to 

80 fm in the ^l^lj calculation with fko = 6.5 MeV. 
dE 



-P 



4 Two-neutron halo in 6 He nucleus 



The 6 He nucleus is studied in the three-body cluster model 6 He = a + n + n. 
The 6 He two-neutron separation energy is 0.976 MeV [41] that is much less 
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than the excitation energy of the a particle lowest excited state and the a 
particle disintegration threshold. As it was already noted, 6 He is a Borromean 
nuclei, i. e. none of the two-body subsystems a + n and n + n has a bound 
state. 

Generally our approach in the 6 He case is very close to the one used in the 
n Li studies discussed above. However, contrary to the n— n Li interaction, 
the n—a interaction is well-known. Therefore it is reasonable to use a more 
microscopically justified approach based on available more realistic n— cluster 
and n—n potentials. In particular, it is reasonable to explicitly allow for the 
spin variables in our model. We couple the spins of two valent neutrons into 
the total spin S — 0, 1 (the a particle has a zero spin), the orbital angular 
momenta l x and l y are coupled to the total orbital angular momentum L, and 
L and S are coupled to the total angular momentum J. Hence the channel 
index T = {K, l x , l y , L, S, J}, and we introduce the following generalization of 
the basis functions (42): 

\KKlJ y (L)(S) : JM) = J2 (LM L SM S \JM) \nKl x l y : LM L ) \SM S ); (53) 

M L , Ms 

the 6 He wave function is given by the general formula (12). In Eq. (53), \SM S ) 
is the spin component of the wave function and \kKI x I v : LMf) is given 
by Eqs. (42)-(43). We note here that allowing for the triplet (S = 1) spin 
states, we enlarge essentially the number of basis functions with any given 
number of oscillator quanta N = 2k + K . Hence for any given N, we have 
the truncated Hamiltonian matrix of a much larger rank, and we are able to 
perform calculations with smaller N than in the n Li case. 

It would be very interesting to perform the 6 He studies based on modern so- 
called realistic NN potentials derived from the meson exchange theory and 
perfectly describing NN scattering data and deuteron observables. However 
the oscillator basis matrices of such potentials are extremely large and cannot 
be handled in calculations of lightest nuclei; the realistic potential matrix ele- 
ments decrease slowly with the number of oscillator quanta iV and truncation 
of the Hamiltonian matrix results in the slow convergence of the results. In 
practical applications usually effective interactions are used in calculations. 
Microscopic ab initio approaches (see, e. g., [42]) involve realistic effective in- 
teractions and effective operators derived from the realistic NN potentials. 
Unfortunately the J-matrix formalism (and other formalisms allowing for the 
continuum spectrum effects) is not developed still for the case of ab initio 
models based on realistic effective NN interactions. In this contribution, we 
are interested in the continuum spectrum effects and S'-matrix pole corrections 
to the binding energy, and so we use phenomenological effective potentials for 
n—n and n—a interactions. A large number of phenomenological effective n—n 
and n—a potentials is available. We perform calculations with various n—n 
and n—a potentials since it is not clear from the very beginning what is the 
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best choice of these potentials; we would like also to compare our results with 
the results of other authors who used different combinations of these interac- 
tions. 

The following effective n—n potentials were employed in our studies. 

The Gaussian potential (38) with Vq — 30.93 MeV and R = 1.82 fm in the 
singlet (S = 0) state and V = 60.9 MeV and R = 1.65 fm in the triplet (S = 1) 
state suggested in Ref. [43], is hereafter refered to as Gs. The potential was 
fitted to the s wave nucleon-nucleon phase shifts only; the interaction between 
the neutrons in the states with relative orbital angular momentum l x > 0, is 
neglected. 

We also make use of a bit different Gaussian potential suggested in Refs. [44,45] 
with V Q = 31.00 MeV and R = 1.8 fm in the singlet state and V = 71.09 MeV 
and R = 1.4984 fm in the triplet state. The singlet component of this potential 
was employed in our n Li studies. We shall refer to this potential as G. 

A more realistic Minnesota n—n potential of Ref. [46] will be refered to as MN. 
This Gaussian potential includes central, spin-orbit and tensor components. 

The lowest single particle 0si/ 2 is occupied in the a particle. There are two 
conventional approaches to the problem of the Pauli forbidden 0si/ 2 state in 
the n + a system. The first approach is to add a phenomenological repulsive 
term to the s wave component of the n— a potential. This phenomenological 
repulsion excludes the Pauli forbidden state in the n+a system and is supposed 
to simulate the Pauli principle effects in more complicated cluster systems. 
This idea was utilized in the SBB n—a potential suggested in Ref. [47]. This 
is an /-dependent Gaussian potential that includes the spin-orbit component. 

Another approach is to use a deep attractive n—a potentials that support the 
Pauli forbidden 0«i/ 2 state in the n + a system. The Pauli forbidden states 
should be excluded in the three-body cluster system. The conventional method 
is to supplement the deep attractive n—a potential by the projecting pseudo 
potential (see, e. g. [48]) 

A|0s 1/2 )(0s 1/2 |. (54) 

If the parameter A is positive and large enough, the projector (54) pushes 
the Pauli forbidden states to very large energies and cleans up the ground 
and low-lying states from the Pauli forbidden admixtures. The eigenfunction 
of the Pauli forbidden state supported by the deep attractive n—a potential 
should be used as |0si/ 2 ) (the so-called eigen-projector [48]); in this case the 
pseudo potential (54) does not affect the description of the scattering data 
provided by the initial deep attractive n—a potential. 

The deep attractive Woods-Saxon n—a potential suggested in Ref. [49] will 
be hereafter refered to as WS. We use the WS potential parameters suggested 
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in Ref. [36] where the radius of the potential was increased in order to fit the 
6 He binding energy. 

We use also deep attractive Majorana splitting potential suggested in Ref. [48] 
(hereafter refered to as MS) and deep attractive /-dependent Gaussian poten- 
tial proposed in Ref. [26] (hereafter refered to as GP). These potentials improve 
the description of scattering data in high (I > 1) partial waves. 

The convergence patterns in the 6 He case are very similar to the ones discussed 
in the n Li case. As a typical example, we present the results obtained in the 
Gs + SBB potential model. The 6 He ground state energy dependence on the 
oscillator parameter fkj for a number of truncation boundary iV values, is 
depicted in Fig. 7. The hio ground state energy dependence is seen to be of 
the same type as the one in the n Li case shown in Fig. 1. The 6 He ground 
state energy and rms radius convergence with N is presented in Fig. 8 for 
two huj values. As in the n Li case (see Fig. 2), there is the staggering of the 
J-matrix binding energy and rms radius with N for frw = 10 MeV that is close 
to the value providing the best convergence of the variational calculation, and 
there is the variational-type iV dependence of these observables for larger Huj 
values. It is interesting that all calculations shown on the right panel of Fig. 7, 
result in nearly the same binding energy at two ~ 15 MeV; this huj value is 
clearly the best convergence value for the J-matrix calculation and it differs 
essentially from the best convergence hu value for the variational calculation 
corresponding to the minima of the curves on the left panel of Fig. 7. 

The results of our calculations of the 6 He two-neutron separation energy and 
rms matter radius, are summarized in Table 2. Since the calculations using 
different potential models were performed with slightly different hu values, we 
list these values in the Table. If for a given potential model the theoretical 
predictions of other authors using other approximations within the three-body 
cluster model are available, they are presented in the corresponding rows of the 
Table, too. The theoretical predictions within the three-body cluster model 
with other potential models and the experimental data, are presented in the 
Table in additional rows. 

The J-matrix approach is seen from the Table to improve the variational 
results for both the binding energy and the rms radius in the case of any 
potential model. The improvement is generally more essential in the case of 
smaller binding energy (it is also seen comparing the 6 He results in Table 2 
with the n Li results in Table 1). The J-matrix approach is also seen to have 
a faster convergence than the variational one, and again the improvement of 
the convergence rate is more essential in the smaller binding energy case. Our 
results are in good correspondence with the available results of other authors 
who used the same potential models. 
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Fig. 7. The 6 He ground state energy in the Gs+SBB potential model vs the oscillator 
basis parameter Huj in the variational (left panel) and J-matrix (right panel) cal- 
culations for different values of the truncation boundary N. The horizontal dashed 
line depicts the convergence limit for the ground state energy. 
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Table 2 

The 6 He two-neutron separation energy and rms matter radius obtained in various potential models and available three-body cluster 
model results of other authors. 



Potential model 


nw, Mev 


Approximation 


E(2n), MeV 
N = 26 N = 28 


( r 2)l/2 > fm 

N = 26 N = 28 


Gs + SBB 


1U.UU 

Re 


variational 
J-matrix 

f. [47] 


0.917 0.919 
0.927 0.923 

"correct asymptotic value" 


2.502 2.510 
2.554 2.539 


Gs + WS 


11.90 

Re 


variational 
J-matrix 

f. [36] 


0.993 0.997 
1.003 1.008 

1.00 


2.394 2.404 
2.444 2.461 

2.44 


G + GP 


11.35 


variational 
J-matrix 


0.875 0.878 
0.889 0.883 


2.386 2.393 
2.451 2.428 


MN + GP 


10.75 


variational 
J-matrix 


0.494 0.509 
0.515 0.516 


2.477 2.484 
2.588 2.547 


MN + MS 


13.00 


variational 
J-matrix 


0.656 0.672 
0.684 0.681 


2.382 2.389 
2.482 2.445 


other potential models 


Ref. [48] : 0.305; Ref. [50] : 1.00; 
Ref. [511: 0.784; Ref. [521: 0.696 


Ref. [50]: 2.50 


experimental data 


Ref. [41]: 0.976 


Ref. [53] 
Ref. [54] 
Ref. [55] 


2.33±0.04 
2.48±0.03 
2.57±0.10 



The structure of the 6 He ground state wave function is presented in Table 3. 
We note here that within the three-body cluster model 6 He = a + n + n, 
the + ground state can be obtained by the coupling of the Jacobi orbital 
momenta l x and l y to the total orbital angular momentum L and subsequent 
coupling of L with the total neutron spin S into the total angular momentum 
J = 0, if only L = S and l x = l y . We arrange the dominant components 
in a different manner. In the Table 3 (a) we show the total weights of all 
components with given S = L and l x = l y , i. e. we sum the contributions 
of the components with different K and N = 2k + K values for the given 
S = L and l x = l y values. In the Table 3 (b) we collect the total weights 
of all components with given S = L and hypermomentum K summing the 
contributions with different l x = l y values. It is also interesting to calculate 
the shell model type component weights which are characterized by the orbital 
momenta l nia and l n2a of individual neutron motion relatively to the a core. 
The unitary transformation from the hyperspherical basis (53) to the shell 
model basis can be found in textbooks (see, e. g., [59,60]). The 6 He + ground 
state wave function can be arranged if only l nia = l n2a . The shell model-type 
component weights are listed in the Table 3 (c). 

It is seen from Table 3 that different potential models result in the wave 
functions of nearly the same structure. For example, the component weights 
in different arrangements obtained with SBB potential simulating the Pauli 
principle by repulsive terms and with WS potential when the Pauli forbidden 
states are explicitly projected out, are nearly the same. Our results for com- 
ponent weights are in good correspondence with the results of other authors 
who used other approaches to the three-body cluster model and other poten- 
tial models. Therefore the uncertainties of the 6 He wave function structure due 
to the uncertainties of the two-body interactions, are very small. The most es- 
sential (however not very large) difference between the dominant component 
weights obtained in different approaches is the difference between the results 
of Ref. [48] where the RSC + SBB* potential model was employed and our 
results and the results of Ref. [47,16] obtained with different potential models. 
We note here that the 6 He binding energy obtained in Ref. [48] is 0.3046 MeV 
only, i. e. it was essentially underestimated. Most probably this is indication 
that the RSC + SBB* potential model is not adequate for the description of 
6 He in the cluster model and hence it is not surprising that the component 
weights of Ref. [48] differ from the ones obtained with other potentials. 

From the naive shell model considerations it follows that two valent neurons 
should occupy p states in the 6 He nucleus. We see from the Table 3 (c) that 
this situation is really utilized in the dominant component of the 6 He wave 
function. The same naive shell model considerations bring us to the K = 2 
component dominance, and again we see from the Table 3 (b) that this is really 
the case. However from these naive considerations it looks strange that the 
dominant component corresponds to the Jacobi orbital momenta l x — l y — 
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Table 3 

Dominant component weights in the 6 He + ground state wave function obtained with various potential models in the J-matrix ap- 
proximation with N = 28 (the respective fko values can be found in Table 2) and the results of other authors in the three-body cluster 
model. 

(a) Total contribution of the components with given S = L and L x = l y (summation over possible K and N values) 



State 








Weight, 


% 






S — L 1% 


— ly 


Gs + SBB 


present work 
Gs + WS G + GP 


MN + MS 


Ref. [47,16] 
Gs + SBB GPT 1 + SBB 


Ref. [48] 
RSC 2 + SBB* 3 







82.1 


83.6 


80.2 


76.6 




83.13 


82.87 


88.908 





2 


2.46 


1.90 


1.64 


2.07 




1.77 


2.31 


1.035 




4 


0.18 


0.14 


0.11 


0.15 




0.02 


0.58 






1 


14.3 


13.4 


17.6 


20.1 




14.54 


13.96 


9.692 


1 


3 


0.85 


0.83 


0.66 


0.89 




0.55 


0.69 


0.366 




5 


0.09 


0.09 


0.06 


0.09 






0.06 





1 Gogny-Pires-de Tourreil n—n potential [56]. 
2 Reid soft-core n—n potential [57]. 

3 SBB n— a potential with modified parameters (see Ref. [48]). 



Table 3 
Prolongation 



(b) Total contribution of the components with given K and S = L (summation over possible L x = l y values) 



State 


Weight, % 


S = L 


K 


Gs + SBB 


present work 
Gs + WS G + GP 


MN + MS 


Ref. [47,16] 
Gs + SBB GPT + SBB 







4.28 


4.06 


3.99 


3.93 


4.41 


4.68 




2 


77.2 


78.9 


75.5 


72.4 


78.93 


78.10 





4 


0.59 


0.19 


0.13 


0.16 


0.53 


0.64 




6 


1.60 


1.61 


1.40 


1.73 


1.16 


1.43 




8 


0.67 


0.51 


0.43 


0.41 


0.39 






2 


14.2 


13.2 


17.5 


20.0 


13.91 


13.41 




4 


0.12 


0.15 


0.05 


0.06 


0.12 


0.15 


1 


6 


0.82 


0.80 


0.63 


0.85 


0.53 


0.67 




8 


0.05 


0.04 


0.04 


0.06 


0.02 





Table 3 
Prolongation 



(c) Total contribution of the components with given shell model-like orbital angular momenta l nia = I 



State 


Weight, % 




Gs + SBB 


present work 
Gs + WS G + GP 


MN + MS 


Ref. [58] 
GPT + WS + V 3 4 


S 


8.27 


7.60 


7.35 


7.26 


7.7 


P 


90.4 


90.9 


91.5 


91.6 


91.0 


d 


0.57 


0.64 


0.43 


0.44 




f 


0.57 


0.54 


0.47 


0.41 




9 


0.07 


0.07 


0.07 


0.05 





4 Three-body nna potential, see Ref. [58]. 



'He (G+GP) 




2 4 6 8 10 12 



E, MeV 

Fig. 9. Reduced El transition probability d in 6 He obtained in the G + GP 
potential model. Vertical lines with cross at the end, dashed and solid lines were 
obtained in the VV, VJ and JJ approximations, respectively, with TvjJ = 11.35 MeV 
and N = 22 for the ground state and N = 23 for the final state. 

[see the Table 3 (a)]. Nevertheless there is, of course, no contradiction between 
the results presented in the Table 3 (a) and the Table 3 (c), they are obtained 
by the summation over different quantum numbers of the same 6 He ground 
state wave functions. 



An example of the reduced (cluster) El transition probability calculations 
in 6 He is shown in Fig. 9 where the results obtained with G + GP potential 
model in the VV, VJ and JJ approximations are presented. In the 6 He case, 
the transition strength is distributed in the VV approximation over a number 

of strong 5-peaks. Therefore the shape of the reduced probability - 

dE 

appears to be more complicated than in the n Li case when the effects of 
continuum are allowed for in the VJ and JJ approximations. As in the n Li 
case, we see that the external asymptotic part of the model space allowed for 
in the JJ approximation only, provides a very significant contribution to the 
electromagnetic transition probabilities. 

The convergence of the reduced El transition probability calculations in the 
J J approximation is illustrated by Fig. 10 where we show the results obtained 
with different truncation boundaries in the ground state calculations N g . s ., 
the truncation boundary in the final state calculations iVf B = N gs +1. The 
dB(El) 

~d~E 

For example, the weakest transition strength is obtained with jV g . s . = 12 and 
Nf.s. = 13 shown by the dotted line in the figure; the strongest transition 
strength is obtained with N gmBm = 14 and iVf s = 15 shown by the dash-dot 
line; the transition strength obtained with iVg. s . = 16 and A f s . = 17 is stronger 
than the one obtained with 7V g . s . = 12 and JVf. s . = 13 but weaker than all the 



staggering of the — — — shape as iV g . s . increases, is clearly seen in the figure. 
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Fig. 10. Convergence of the dt3 }^ cal- 
culations in 6 He with the truncation pa- 
rameter iVg. s . used in the ground state 
calculations (the final state truncation 
parameter 7Vf_ s _ = N g . s . + 1). The re- 
sults are obtained in the JJ approxima- 
tion with the G + GP potential model 
and Huj = 11.35 MeV. Arrows show the 
changes of d when N g . s . is increased 
by 4. 



Fig. 11. Reduced El transition proba- 
bility — -j-gp- in 6 He obtained with vari- 
ous potential models in the JJ approxi- 
mation with the ground state truncation 
parameter N g _ s _ = 22 and the final state 
truncation parameter iVf. s . = 23 (the 
corresponding hio values can be found 
in Table 2) in comparison with the re- 
sults of Funada et al. [61] and Danilin 
et al. [62]. 



rest results presented in the figure, etc. However selecting the results obtained 
with even N g _ s j2 values only, we see that the El transition strength increases 
monotonically with N g . s .', selecting the results obtained with odd N g _ s _/2 values 
only, we see the monotonic decrease of the El transition strength with N g . s .- 

The comparison of the reduced El transition probability results obtained with 
different potential models and calculations of other authors within the three- 
body cluster model, is shown in Fig. 11. 



5 Phase equivalent transformation with continuous parameters and 
three-body cluster system 



The results presented in Table 2 show that we obtain a very good approxima- 
tion for the 6 He two-neutron separation energy. However the exact value of the 
6 He binding energy is not reproduced with any potential model employed. We 
suggest to use a phase equivalent transformation of the two-body interaction 
to improve the description of the 6 He binding. If the phase equivalent trans- 
formation depends on a continuous parameter(s) than varying the parameter 
we can fit the 6 He binding energy to the phenomenological value. We suppose 
this approach to be interesting for various few-body applications. 
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Various phase equivalent transformations have been discussed in literature. 
The local phase equivalent transformations (transforming a local potential 
into another local potential phase equivalent to the initial one) is well-known 
(see, e. g., Ref. [63]). In the inverse scattering theory, this transformation 
gives rise to the ambiguity of the potential restored from the scattering data. 
Recently this transformation was extended on the case of the multichannel 
scattering [64]. However the local phase equivalent transformation can be ap- 
plied only to a two-body system that have at least one bound states and the 
number of continuous parameters of the transformation is equal to the number 
of the bound states in the system. Therefore this transformation cannot be 
applied to the n + a system that has no bound states. 

Recently the so-called supersymmetry phase equivalent transformation (see, 
e. g., review [65]) become very popular. Using this transformation one can 
transform a potential with the Pauli forbidden state into exactly phase equiv- 
alent potential with additional repulsion simulating the Pauli effects. The ef- 
fect of this transformation on the properties of three-body cluster systems was 
examined in a number of recent papers (see, e. g., Ref. [66]). The supersym- 
metry transformation was shown [64] to be a particular case of the local phase 
equivalent transformation. However the supersymmetry transformation does 
not have parameters and cannot be used for our purposes. 

A phase equivalent transformation based on unitary transformation of the 
Hamiltonian, was suggested in Ref. [67]. This transformation has continuous 
parameters. This transformation have been used in many-body applications 
in Ref. [68]. However the authors of Ref. [68] used few particular cases of the 
transformation corresponding to particular parameter values and did not try 
to vary the parameters continuously. 

We have developed recently [26] a phase equivalent transformation using the 
J-matrix formalism. Generally, our transformation is a particular case of the 
phase equivalent transformation of Ref. [67]. However we suppose that our 
phase equivalent transformation is general enough and very convenient for the 
use in various many-body applications utilizing any L 2 basis. We applied this 
transformation to the NN interaction and studied the effect of the transfor- 
mation on the properties 3 H and 4 He nuclei in Ref. [69]. In what follows, we 
discuss briefly the transformation and its application to the 6 He nucleus within 
the three-body cluster model. 

We suppose that a two-body system is described by the Schrodinger equation 

HV(r) = EV(r), (55) 

where the Hamiltonian H = T + V, T and V are kinetic and potential energy 
operators. Introducing a complete basis {0 K (?")}, k — 0, 1, 2, ... of L 2 func- 
tions <f> K (r), we can expand the solutions of Eq. (55) in infinite series of basis 
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functions \k) = K (r): 



*(r) = E^I«> 



(56) 



k=0 



The Schrodinger equation (55) takes the form of an infinite dimensional alge- 
braic problem 

oo 

J2(k\H\k')C k ,=EC k , (57) 

k'=0 

where (k\H\k') are the matrix elements of the infinite dimensional Hamiltonian 
matrix [H\. 

Now we define a new matrix 

[H'\ = [[/+] [H] [U] (58) 
with the help of the unitary matrix [U] which is supposed to be of the form 

/ 



[U] = [U ] © [/] 











[I] 



(59) 



where [J] is the infinite dimensional unit matrix and [Uq] is N x N unitary ma- 
trix. A new Hamiltonian H' is defined through its matrix [H'\. It is supposed 
that [H'\ is the matrix of the Hamiltonian H' in the original basis {4> K { r )}- 

Clearly the spectra of Hamiltonians H and H' are identical. The difference 
between the eigenfunctions \l/'(r) of the Hamiltonian H' and the eigenfunc- 
tions ^(r) of the Hamiltonian H corresponding to the same energy E, is a 
superposition of a finite number of square integrable functions: 

N-l 

*'(r) = W(r)+ ]T AC K K (r). (60) 

ft=0 

The superposition of a finite number of L 2 functions cannot affect the asymp- 
totics of scattering wave functions. Since the scattering phase shifts and the 
S'-matrix are defined through the asymptotic behavior of the wave functions, 
the phase shifts associated with the functions \&(r) and ^'(r) are identical. In 
other words, the Hamiltonians H and H' are phase equivalent. 

Supposing that the Hamiltonian H' = T + Vpet, we introduce the potential 

V PET = V + AV (61) 

phase equivalent to the original potential V. The potential Vpet is defined 
through its matrix 



[Vpet] = [V] + [AV] 



(62) 
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in the basis {0 K (r)}, where 



[H] [U] - [H] . (63) 

The net result of the above considerations is very simple. We introduce any 
complete L 2 basis {(f> n {^)} and calculate the Hamiltonian matrix in this basis. 
Next we introduce an arbitrary unitary transformation of the type (59) that 
affects a finite number of the basis functions. With the help of this transforma- 
tion, we calculate the additional potential AV using Eq. (63) and obtain the 
phase equivalent interaction by means of Eq. (61). The additional potential 
AV is non-local, and hence the phase equivalent potential Vpet is non-local, 
too. Therefore the suggested transformation is a non-local phase equivalent 
transformation: it brings us to a non-local interaction phase equivalent to the 
original one. 

The suggested non-local phase equivalent transformation can be easily imple- 
mented in many-body calculations utilizing any L 2 basis. One just needs to 
add the two-body kinetic energy matrix to the two-body interaction matrix, 
unitary transform the obtained matrix and subtract the kinetic energy ma- 
trix from the result to obtain the matrix of the phase equivalent two-body 
interaction. 

The non-local phase equivalent transformation can be easily explained and 
understood within the J-matrix formalism. In the J-matrix formalism, the S- 
matrix and the phase shifts are governed by the matrix elements 
(KrT\^\Kr' + 1, T') which are defined through eigenvalues E\ and the last 
component (/tpTlA) of the eigenvectors of the truncated Hamiltonian matrix 
[see Eqs. (30)-(31)]. If the truncated Hamiltonian matrix is larger than the 
non-trivial submatrix [Uq] of the unitary matrix [U], than neither E\ nor 
the last component (/c r r|A) of the eigenvectors are affected by the unitary 
transformation (59). Therefore neither the S'-matrix nor the phase shifts are 
affected by the transformation. Nevertheless the wave function is seen from 
Eq. (60) to be subjected to changes by the transformation; in other words, the 
off-shell properties of the original and the transformed potentials are different. 
The non-local phase equivalent transformation generates the ambiguity of the 
interaction obtained by means of the J-matrix version of the inverse scattering 
theory [70]. 

The off-shell properties of the two-body interactions play an important role 
in the formation of the properties of a three-body (and generally many-body) 
system. Therefore it is interesting to investigate the effect of the non-local 
phase equivalent transformation on the properties of the 6 He nucleus within 
the three-body cluster model. We note that the non-local phase equivalent 
transformation preserves the energies of the two-body bound states. However 
the properties of these states are not preserved, for example, the rms radius 



[AV] 
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of the two-body bound system may be affected by the transformation [71]. 
Hence it is somewhat dangerous to apply the transformation to the NN in- 
teraction since it may destroy the description of the deuteron properties. We 
have already noted that we employ effective n—n interactions that are de- 
signed for the effective description of many-body nuclear systems and that are 
not supposed to be careful in the description of the dinucleon. Nevertheless 
we suppose that it is more natural to apply the transformation to the n— a 
interaction since the n + a system does not have a bound state and all the 
information about this interaction is extracted from the scattering data only. 
It is clear from the Table 3 (c) that the 6 He properties are dominated by the 
p wave component of the n— a interaction. Therefore we apply the non-local 
phase equivalent transformation to the p wave component of the n— a potential 
only. 

The non-local phase equivalent transformation for the n— a potential p wave 
component is constructed in the oscillator basis with tko values used in three- 
body calculations with the respective potential model (see Table 2). The sim- 
plest non-local transformation involves a 2 x 2 unitary matrix [Uq\. Any 2x2 
unitary matrix is known to be a rotation matrix with a single continuous 
parameter 7: 



(64) 



cos 7 — sin 7 
+ sin 7 cos 7 

A more complicated non-local transformation involves a 3 x 3 rotation matrix 
[Uq] with two continuous parameters 7 and j3: 



cos (3 cos 7 — cos P sin 7 sin f3 

sin 7 cos 7 

— sin p cos 7 sin ft sin 7 cos f3 



(65) 



Clearly the transformation with the matrix (65) is equivalent to the transfor- 
mation with the matrix (64) if the Euler angle /3 — 0. 

We present in Fig. 12 the results of variational calculations of the 6 He ground 
and the first excited + states with A = 20 and of the lowest 1~ state with 
A = 21 with various n—n potentials and n— a potentials obtained by applying 
the phase equivalent transformation (64) to various original potentials. The 
7 dependence of + and 1~ state energies is very interesting. It is seen from 
Fig. 12 that the simplest non-local phase equivalent transformation of p wave 
component of the n— a potential only, can completely change the spectrum of 
the three-body cluster system. The 7 dependence of the energies is seen to be 
very similar for all potential models under consideration. The variations of the 
energies of the excited + and l - states are seen to be smaller than the ground 
state energy variations. There is the + states level crossing at 7 —150° for 
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Fig. 12. The 6 He ground state energy (solid line) and the energies of the lowest 
1 _ (dashed line) and + (dotted line) excited states obtained in the variational 
approach with phase equivalent n—a potentials vs the transformation parameter 7 
of the one-parameter transformation (64). The truncation parameter N = 20 in the 
+ state calculations and N = 21 in the 1~ state calculations. 



all potential model considered. 



The ground state energy 7 dependence passes through a minimum in the vicin- 
ity of 7 = 5— 7°. The minimum corresponds to the increase of the binding by 
approximately 0.04—0.07 MeV or 7—12% depending on the potential model, 
that improves the results obtained with the original potentials (7 = 0) given 
in Table 2. The ground states may be additionally shifted down using the two- 
parameter transformation (65). Supposing 7 = 7° and varying the parameter 
P we obtain variational results shown in Fig. 13. The ground state (3 depen- 
dence has a minimum at /3 ~ —2.5° for all potential models. This minimum 
corresponds to the additional binding of 0.019 -h 0.023 MeV (about 2-4%). 

The variational results presented in Figs. 12 and 13 are interesting for under- 
standing the general trends of variation of the ground state energy and 6 He 
spectrum when the parameters of the phase equivalent transformations (64)- 
(65) are varied in a wide range of values. The corrections to the 6 He binding 
due to the effect of the phase equivalent transformation are better seen in 
Fig. 14 where we present in a larger scale the J-matrix results for the ground 



34 







6 He 






7 = 7° 




3 












MeV 


2 


■/ 












1 


•' / 


\l 

It 
U 
11 

\l 















U 
U 


i 

t 










n 
\i 


el 
ci 








1 














-180 -120 


-60 60 


120 180 








|8, deg. 





J— matrix 



> 




Fig. 13. The 6 He ground state en- 
ergy vs the transformation parameter (3 
of the phase equivalent transformation 
(65) with 7 = 7°. Solid, dashed and dot- 
ted lines are calculations in the varia- 
tional approximation with G + GPpet, 
MN + GPpet and MN + MS PE t poten- 
tial models, respectively, with Hlo values 
given in Table 2 and N = 20. 



Fig. 14. The 6 He ground state en- 
ergy obtained in the J-matrix approach 
with phase equivalent n—a potentials vs 
the transformation parameter 7 of the 
one-parameter transformation (64). The 
truncation parameter N = 20, the hw 
values for each potential model can be 
found in Table 2. The straight dashed 
line is the experimental ground state en- 
ergy. 



state energies. The 6 He binding energy is seen to be very close to the empirical 
value in the G + GPpet potential model when 7^7°. However the results 
presented in Fig. 14 were obtained with not very large value of the truncation 
parameter TV = 20. With N = 28 and 7 = 7.5° we obtain in this potential 
model an excellent description of the 6 He binding energy E b = 0.952 MeV and 
rms radius (r 2 ) 1 / 2 = 2.37 fm. Using the two-parameter transformation (65) 
with 7 = 7.5° and (5 = -3°, we obtain E b = 0.973 MeV and (r 2 ) 1/2 = 2.36 fm, 
i. e. the phenomenological value E b = 0.976 MeV is reproduce nearly exactly. 

The correlations between the 6 He rms radius (r 2 ) 1 / 2 and the square root of 

1/2 

the binding energy E b , are depicted in Fig. 15. We present the J-matrix 
results obtained with various potential models when the parameter 7 of the 
phase equivalent transformation (64) varies on the interval where 6 He appears 
to be bound in the variational approximation in the given potential model. 
The correlations are seen to be very interesting: there are two very different 
(r 2 ) 1 / 2 values that are in correspondence with the same binding energy. For 
7 < —6°, the rms radius (r 2 ) 1 / 2 decreases linearly with E\j 2 . 

The effect of the phase equivalent transformation (64) on the reduced El 

transition probability — - in 6 He, is illustrated by Fig. 16. The effect is 

seen to be essential. The naive expectation is that the soft dipole mode will be 
more enhanced if the binding energy is smaller. However the strength of the El 



35 



6 He 




2.2 >— ' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
0.0 0.2 0.4 0.6 0.8 1.0 

E{ /z , MeV 1/2 

Fig. 15. Correlations between E^J 2 and rms radius in 6 He obtained by variation 
of the parameter 7 of the phase equivalent transformation (64) in various potential 
models ( J-matrix approximation with N = 20, the corresponding hio values can be 
found in Table 2). A straight solid line is added to visualize the linear correlation 
between (r 2 ) 1 / 2 and E^ 2 on a part of the trajectory. 



. , . . . r , . r dB(El) , 

transitions m the vicinity 01 the maximum 01 — — — — does not demonstrate 

d E 

so simple dependence on the binding energy E b (the corresponding E b values 
are listed in the figure). 
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Fig. 16. Reduced El transition probability dB }^ in 6 He obtained in the JJ ap- 



proximation with the ground state truncation parameter iV g . s . = 22 and the fi- 
nal state truncation parameter 7Vf. s . = 23 with the MN + GPpet (left panel) and 
MN + MSpet (right panel) potential models (the corresponding Huj values can be 
found in Table 2) and different values of the parameter 7 of the phase equivalent 
transformation (64) (the corresponding binding energies are listed in the figure). 
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6 Conclusions 

The hyperspherical J-matrix formalism makes it possible to study not only 
the few-body disintegration of the system but also to calculate the ^-matrix 
poles in the A-body system. This approach appears to be a very powerful tool 
for calculations of the bound state properties. As a result, we obtain a unified 
theory that is capable to investigate in a unique approach both the discrete 
and continuum spectra of A-bodj systems. 

The suggested J-matrix motivated non-local phase equivalent transformation 
may be used to fit binding energies of many-body systems in the case when all 
information about two-body interaction is extracted from the scattering data 
only. The transformation can be easily utilized in the studies of many-body 
systems with any L 2 basis. 

This work was supported in part by the State Program "Russian Universities" 
and by the Russian Foundation of Basic Research, Grant No 02-02-17316. 



References 



[1] E. J. Heller and H. A. Yamani, Phys. Rev. A 9, 1201 (1974); 9, 1209 (1974). 

[2] H. A. Yamani and L. Fishman, J. Math. Phys. 16, 410 (1975). 

[3] J. T. Broad and W. P. Reinhardt, Phys. Rev. A 14, 2159 (1976); J. Phys. B 9, 
1491 (1976). 

[4] A. M. Shirokov, Yu. F. Smirnov, and L. Ya. Stotland, in Proc. Xllth Europ. 
Conf. on Few-Body Phys., Uzhgorod, USSR, 1990 (Ed. V. I. Lengyel and 
M. I. Haysak), 173 (Uzhgorod, 1990). 

[5] Yu. F. Smirnov, L. Ya. Stotland, and A. M. Shirokov, Izv. Akad. Nauk SSSR, 
Ser. Fiz. 54, No 5, 897 (1990) [Bull. Acad. Sci. USSR, Phys. Ser. 54, No 5, 81 
(1990)]. 

[6] D. A. Konovalov and I. A. McCarthy, J. Phys. B27, L407 (1994); J. Phys. B27, 
L741 (1994). 

[7] G. F. Filippov and I. P. Okhrimenko, Yad. Fiz. 32, 932 (1980) [Sov. J. Nucl. 
Phys. 32, 480 (1980)]; G. F. Filippov, Yad. Fiz. 33, 928 (1981) [Sov. J. Nucl. 
Phys. 33, 488 (1981)]. 

[8] Yu. F. Smirnov and Yu. I. Nechaev, Kinam 4, 445 (1982); Yu. I. Nechaev and 
Yu. F. Smirnov, Yad. Fiz. 35, 1385 (1982) [Sov. J. Nucl. Phys. 35, 808 (1982)]. 



37 



[9] G. F. Filippov, V. S. Vasilevski, and L. L. Chopovski, Fiz. Elem. Chastits At. 
Yadra 15, 1338 (1984); 16, 349 (1985) [Sov. J. Part. Nucl. 15, 600 (1984); 16, 
153 (1985)]. 

[10] G. F. Filippov, Rivista Nuovo Cim. 12, 1 (1989). 

[11] V. A. Knyr, A. I. Mazur, and Yu. F. Smirnov, Yad. Fiz. 52, 754 (1990) [Sov. 
J. Nucl. Phys. 52, 483 (1990)]. 

[12] V. A. Knyr, A. I. Mazur, and Yu. F. Smirnov, Yad. Fiz. 54, 1518 (1991) [Sov. 
J. Nucl. Phys. 54, 927 (1991)]. 

[13] S. P. Merkuriev and L. D. Faddeev, Quantum scattering theory for systems of 
few bodies (Nauka Publishes, Moscow, 1985). 

[14] R. I. Jibuti, Elem. Part. Atom. Nucl. 14, 741 (1983). 

[15] R. I. Jibuti and N. B. Krupennikova, Hyper spherical harmonics method in 
quantum mechanics of few bodies (Metsniereba, Tbilisi, 1984). 

[16] M. V. Zhukov, B. V. Danilin, D. V. Fedorov, J. M. Bang, I. J. Thompson, and 
J. S. Vaagen, Phys. Rep. 231, 151 (1993). 

[17] A. M. Shirokov, in Proc. Int. Conf. Perspectives of Nuclear Physics in the 
Late Nineties, Hanoi, Vietnam, 1994 (Ed. Nguyen Dinh Dang, Da Hsuan Feng, 
Nguyen Van Giai, Nguyen Dinh Tu), 295 (Singapore, World Scientific, 1995); 
|nucl-th/9407019l (1994). 

[18] J. M. Bang, B. V. Danilin, V. D. Efros, J. S. Vaagen, M. V. Zhukov, and 
I. J. Thompson, Phys. Rep. 264, 27 (1996). 

[19] E. Nielsen, D. V. Fedorov, A. S. Jensen, and E. Garrido, Phys. Rep. 347, 373 
(2001). 

[20] Yu. F. Smirnov and A. M. Shirokov, Preprint ITF-88-47R (Kiev, 1988); A. M. 
Shirokov, Yu. F. Smirnov, and S. A. Zaytsev, in Modern Problems in Quantum 
Theory (Ed. V. I. Savrin and O. A. Khrustalev), (Moscow, 1998), 184; Teoret. 
Mat. Fiz. 117, 227 (1998) [Theor. Math. Phys. 117, 1291 (1998)]. 

[21] T. Ya. Mikhelashvili, Yu. F. Smirnov, and A. M. Shirokov, Yad. Fiz. 48, 969 
(1988) [Sov. J. Nucl. Phys. 48, 617 (1988)]; J. Phys. G 16, 1241 (1990). 

[22] D. E. Lanskoy, Yu. A. Lurie, and A. M. Shirokov, Z. Phys. A 357, 95 (1997). 

[23] Yu. A. Lurie, Yu. F. Smirnov, and A. M. Shirokov, Izv. Ros. Akad. Nauk, Ser. 
Fiz. 57, 193 (1993) [Bull. Rus. Acad. Sci., Phys. Ser. 57, 943 (1993)]. 

[24] A. M. Shirokov, Yu. A. Lurie, and Yu. F. Smirnov, Halo'94 NORDITA /NorFa 
Study Weekend on Halo Nuclei, Copenhagen, 1994 (Ed. J. M. Bang and 
J. S. Vaagen). 

[25] Yu. F. Smirnov, Yu. A. Lurie, and A. M. Shirokov, Revista Mex. Fis. 40, Supl. 1, 
63 (1994). 



38 



[26] Yu. A. Lurie and A. M. Shirokov, Izv. Ros. Akad. Nauk, Ser. Fiz. 61, 2121 
(1997) [Bull. Rus. Acad. Sci., Phys. Ser. 61, 1665 (1997)]. 

[27] I. P. Okhrimenko, Few-body Syst. 2, 169 (1987). 

[28] Yu. F. Smirnov and K. V. Shitikova, Fis. Elem. Chast. At. Yad. 8, 847 (1977) 
[Sov. J. Part. Nucl. 8, 344 (1977)]. 

[29] J. M. Bang, A. I. Mazur, A. M. Shirokov, Yu. F. Smirnov, and S. A. Zaytsev, 
Ann. Phys. (NY) 280, 299 (2000). 

[30] L. Johansen, A. S. Jensen, and P. G. Hansen, Phys. Lett. B 244, 356 (1990). 

[31] W. Benenson, Nucl. Phys. A 588, 11c (1995). 

[32] M. Abramowitz and I. A. Stegun (eds.), Handbook on Mathematical Functions 
(Dover, New York, 1972). 

[33] M. V. Zhukov, B. V. Danilin, D. V. Fedorov, J. S. Vaagen, F. A. Gareev, and 
J. Bang. Phys. Lett. B 265, 19 (1991). 

[34] F. Ajzenberg-Selov, Nucl. Phys. A 490, 1 (1988). 

[35] I. Tanihata, T. Kobayashi, O. Yamakawa, S. Shimoura, K. Ekuni, K. Sugimoto, 
N. Takahashi, T. Shimoda, and H. Sato, Phys. Lett. B 206, 592 (1988). 

[36] J. S. Al-Khalili, J. A. Tostevin, and I. J. Thompson, Phys. Rev. C 54, 1843 
(1996). 

[37] J. S. Al-Khalili and J. A. Tostevin, Phys. Rev. Lett. 76, 3903 (1996). 

[38] H. Sagawa and M. Houma, Phys. Lett. B 251, 17 (1990); H. Sagawa, 
N. Takigawa, and N. van Giai, Nucl. Phys. A 543, 575 (1992). 

[39] B. V. Danilin, M. V. Zhukov, J. S. Vaagen, I J Thompson, and J. M. Bang, 
Surrey University Preprint CNP92/17 (1992). 

[40] D. Sackett D. et al., Phys. Rev. C 48, 118 (1993); K. Ieki K. et al., Phys. Rev. 
Lett. 70, 730 (1993); A. Galonsky A. et al., Michigan State University Preprint 
MSUCL-899 (1993). 

[41] A. H. Wapstra, G. Audi, and R. Hoekstra, At. Data Nucl. Data Tables 39, 281 
(1988). 

[42] P. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. Lett. 84, 5728 (2000); 
Phys. Rev. C 62, 054311 (2000). 

[43] G. E. Brown and A. D. Jackson, The Nucleon-Nucleon Interaction (North- 
Holland, Amsterdam, 1976). 

[44] B. V. Danilin, M. V. Zhukov, J. S. Vaagen, and J. M. Bang, Phys. Lett. B 302, 
129 (1993). 

[45] L. S. Ferreira, E. Maglione, J. M. Bang, I. J. Thompson, B. V. Danilin, 
M. V. Zhukov, and J. S. Vaagen, Phys. Lett. B 316, 23 (1993). 



39 



[46] T. Kaneko, M. LeMere, and L. C. Tang, University of Minnesota Preprint TPI- 
MINN-91/32-T (1991). 

[47] B. V. Danilin, M. V. Zhukov, A. A. Korsheninnikov, and L. V. Chulkov, Yad. 
Fiz. 53, 71 (1991) [Sov. J. Nucl. Phys. 53, 45 (1991)]. 

[48] V. I. Kukulin, V. N. Pomerantsev, Kh. D. Razikov, V. T. Voronchev, and 
G. G. Ryzhikh, Nucl. Phys. A 586, 151 (1995). 

[49] J. Bang and C. Gignoux, Nucl. Phys. A 313, 119 (1979). 

[50] E. Garrido, D. V. Fedorov, and A. S. Jensen, Nucl. Phys. A 617, 153 (1997). 

[51] K. Kato, S. Aoyama, S. Mukai, and I. Ikeda, Nucl. Phys. A 588, 29c (1995); 
S. Aoyama, S. Mukai, K. Kato, and I. Ikeda, Progr. Theor. Phys. 93, 99 (1995). 

[52] E. Hiyama, M. Kamimura, Nucl. Phys. A 588, 35c (1995). 

[53] I. Tanihata, D. Hirata, T. Kobayashi, S. Shimoura, K. Sugimoto, and H. Toki, 
Phys. Lett. B 289, 261 (1992). 

[54] I. Tanihata et al., RIKEN Preprint AF-NP-60 (1987). 

[55] L. V. Chulkov, B. V. Danilin, V. D. Efros, A. A. Korsheninnikov, and 
M. V. Zhukov, Europhys. Lett. 8, 245 (1989). 

[56] D. Gogny, P. Pires, and R. de Tourreil, Phys. Lett. B 32, 591 (1970). 

[57] R. V. Reid, Jr., Ann. Phys. (NY) 50, 411 (1968). 

[58] I. J. Thompson, B. V. Danilin, V. D. Efros, J. S. Vaagen, J. M. Bang, and 
M. V. Zhukov, Phys. Rev. C 61, 024318 (2000). 

[59] O. F. Nemets, V. G. Neudatchin, A. E. Rudchik, Yu. F. Smirnov, and 
Yu. M. Tchuvil'sky, Nucleon associations in atomic nuclei and multi-nucleon 
transfer reactions (Naukova Dumka Publishers, Kiev, 1988). 

[60] M. Moshinsky and Yu. F. Smirnov, The Harmonic Oscillator in Modern Physics 
(Harwood Academic Publishers, Amsterdam, 1996). 

[61] S. Funada, H. Kameyama, and Y. Sakuragi, Nucl. Phys. A 575, 93 (1994). 

[62] B. V. Danilin, T. Rogde, S. N. Ershov, H. Heiberg-Andersen, J. S. Vaagen, 
I. J. Thompson, and M. V. Zhukov, Phys. Rev. C 55, R577 (1997). 

[63] R. G. Newton, Scattering theory of waves and particles, 2nd. ed. (Springer- 
Verlag, New York, 1982). 

[64] A. M. Shirokov and V. N. Sidorenko, Yad. Fiz. 63, 2085 (2000) [Phys. At. Nucl. 
63 1993 (2000)]. 

[65] F. Cooper, A. Khare, and U. Sukhatme, Phys. Rep. 251, 267 (1995). 

[66] E. Garrido, D. V. Fedorov, and A. S. Jensen, Nucl. Phys. A 650, 247 (1999). 



40 



[67] F. Coester, S. Cohen, B. Day, and C. W. Vincent, Phys. Rev. C 1, 769 (1970); 
M. I. Haftel and F. Tabakin, Phys. Rev. C 3, 921 (1971). 

[68] H. C. Pradhan, P. U. Sauer, and J. P. Vary, Phys. Rev. C 6, 407 (1972). 

[69] A. M. Shirokov, J. P. Vary, T. A. Weber, and A. I. Mazur, Izv. Ros. Akad. 
Nauk, Ser. Fiz. [Bull. Rus. Acad. Sci., Phys. Ser.] 65, 1572 (2001). 

[70] S. A. Zaytsev, Teoret. Mat. Fiz. 115, 263 (1998) [Theor. Math. Phys. 115, 575 
(1998)]; 121, 424 (1999) [Theor. Math. Phys. 121, 1617 (1999)]. 

[71] W. N. Polyzou, Phys. Rev. C 58, 91 (1998). 



41 



